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The correlated motion of electrons in the presence of strong orbital fluctuations 
and correlations is investigated with respect to magnetic couplings and excitations 
in an orbitally degenerate ferromagnet within the framework of a non-perturbative 
Goldstone-mode-preserving approach based on a systematic inverse-degeneracy ex- 
pansion scheme. Introduction of the orbital degree of freedom results in a class 
of diagrams representing spin-orbital coupling which become particularly important 
near the orbital ordering instability. Low-energy staggered orbital fluctuation modes, 
particularly with momentum near (7r/2, 7r/2, 0) (corresponding to period 4a orbital 
correlations as in CE phase of manganites involving staggered arrangement of nomi- 
nally Mn^''"/Mn^+ ions, and staggered ordering of occupied Sx^ — r^/3y^ — orbitals 
on alternating Mn^^ sites), are shown to generically yield strong intrinsically non- 
Heisenberg (1 — cosg)^ magnon self energy correction, resulting in no spin stiffness 
reduction, but strongly suppressed zone-boundary magnon energies in the F-X direc- 
tion. The zone-boundary magnon softening is found to be strongly enhanced with 
increasing hole doping and for narrow-band materials, which provides insight into 
the origin of zone-boundary anomalies observed in ferromagnetic manganites. 
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I. INTRODUCTION 

The orbital degree of freedom of the electron has attracted considerable attention in 
recent years due to the rich variety of electronic, magnetic, and transport properties ex- 
hibited by orbitally degenerate systems such as the ferromagnetic manganites, which have 
highlighted the interplay between spin and orbital degrees of freedom in these correlated elec- 
tron systems. ^'^ Orbital fluctuations, correlations, and orderings have been observed in Ra- 
man spectroscopic studies^ of orbiton modes in LaMnOs, polarization-contrast-microscopy 
studies^ of Lao.5Sri.5Mn04, magnetic susceptibility and inelastic neutron scattering studies^ 
of La4Ru20io, and resonant inelastic soft X-ray scattering studies^ of YTiOs and LaTiOa. 
A new detection method for orbital structures and ordering based on spectroscopic imag- 
ing scanning tunneling microscopy is of strong current interest^ in orbitally active metallic 
systems such as strontium ruthenates and iron pnictide superconductors. 

The role of orbital fluctuations on magnetic couplings and excitations is of strong 
current interest in view of the several zone-boundary anomalies observed in spin-wave 
excitation measurements in the metallic ferromagnetic phase of colossal magnetoresis- 
tive (CMR) manganites.-"— The presence of short-range dynamical orbital fluctuations 
has been suggested in neutron scattering studies of ferromagnetic metallic manganite 
Lai_x(Cai_ySry)xMn03. — These observations are of the crucial importance for a quantitative 
understanding of the carrier-induced spin-spin interactions, magnon excitations, and magnon 
damping, and have highlighted possible limitations of existing theoretical approaches. 

For example, the observed magnon dispersion in the F-X direction shows significant 
softening near the zone boundary, indicating non-Heisenberg behaviour usually modeled by 
including a fourth neighbour interaction term J4, and highlighting the limitation of the 
double-exchange model. Similarly, the prediction of magnon-phonon coupling as the origin 
of magnon damping^ and of disorder as the origin of zone-boundary anomalous softening^^ 
have been questioned in recent experiments.—"— Furthermore, the dramatic difference in the 
sensitivity of long-wavelength and zone-boundary magnon modes on the density of mobile 
charge carriers has emerged as one of the most puzzling feature. Observed for a finite range 
of carrier concentrations, while the spin stiffness remains almost constant, the anomalous 
softening and broadening of the zone-boundary modes show substantial enhancement with 
increasing hole concentration.— 
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Theoretically, the role of orbital-lattice fluctuations and correlations on magnetic cou- 
plings and excitations has been investigated within an orbitally degenerate double exchange 
model with an inter-orbital interaction V and the Jahn- Teller coupling.— Based on a strong 
coupling expansion, this approach is restricted to the strong coupling limit V ^ t. The 
flnal calculations for the magnon self energy, carried out in terms of a phenomenological 
parameter, show signiflcant zone-boundary magnon softening only for ferromagnetic orbital 
correlations, and extremely close to the orbital ordering instability. 

If orbital fluctuations have signature effects on magnetic excitations in a ferromagnet 
with orbital degree of freedom, they can be probed indirectly through neutron scattering 
studies. A detailed investigation of the orbital fluctuation magnon self energy is therefore 
of strong current interest, especially with respect to dependence on inter-orbital interaction 
strength, band fllling, and different orbital fluctuation modes. In this paper we will present 
a theory for spin-orbital coupling and magnon self energy, and examine how the correlated 
motion of electrons in the presence of strong orbital correlations near the orbital ordering 
instability influences magnetic couplings and excitations. 

We will employ a diagrammatic approach which allows interpolation in the full range 
of interaction strength from weak to strong coupling. In this approach correlation effects 
in the form of self energy and vertex corrections are incorporated systematically so that 
the Goldstone mode is preserved order by order. Based on a systematic inverse- degeneracy 
expansion scheme,— the approach has been used recently to study spin-charge coupling 
effects, which give rise to strong magnon energy softening, damping, and non-Heisenberg 
behaviour in metallic ferromagnets.— 

The present work will also extend the recent investigation into role of orbital degeneracy 
and Hund's coupling on magnetic couplings and excitations in a band ferromagnet.™ Or- 
bital degeneracy and Hund's coupling were shown to enhance ferromagnetism by strongly 
suppressing correlation-induced quantum corrections to spin stiffness and magnon energies. 
An effective quantum parameter was obtained for determining the magnitude of quantum 
corrections, and the theory was applied to calculate the spin stiffness for a realistic multi- 
orbital system such as iron. We will show here that the spin stiffness remains essentially 
unaffected by the interaction V due to the non-Heisenberg (1 — cosg)^ behaviour of the 
magnon self energy resulting from orbital fluctuations and correlations. 

In manganites, an important role is also played by the cooperative Jahn- Teller distortion 



4 



of ions which hfts the two-fold degeneracy of eg electronic levels of Mn due to a combina- 
tion of orbital geometry and electrostatic repulsion, leading to staggered orbital correlations. 
This is qualitatively similar to the local orbital moment and staggered orbital correlations 
introduced by the inter-orbital density interaction VniaUif^ which relatively pushes up the 
/3 orbital energy if the a orbital density {uia) is more than average, thus self consistently 
lifting the orbital degeneracy. Therefore, orbital correlations and fluctuations due to dynam- 
ical Jahn-Teller distortion can be qualitatively treated in terms of an effective inter-orbital 
interaction. 

The structure of the paper is as follows. Starting with a degenerate two-orbital Hubbard 
model including an inter-orbital interaction V, the first order quantum correction diagrams 
for the irreducible particle-hole propagator are obtained in Section II. As basic ingredients 
in the diagrammatics, spin and orbital fluctuations are briefly discussed in section III. The 
magnon self energy contributions due to orbital fluctuations and spin-orbital coupling are 
obtained in sections IV and V. The interplay between magnetic and charge contributions to 
the spin-orbital interaction vertex is discussed in section VI, and orbital fluctuations near 
(7r/2, 7r/2, 0) are shown to yield strong zone-boundary magnon softening. Extension to finite 
Hund's coupling J and the ferromagnetic Kondo lattice model are discussed in sections VII 
and VIII, and conclusions are presented in Section IX. 

II. TWO-ORBITAL HUBBARD MODEL 

We will consider a two-orbital Hubbard model 

(ij)<T i 
+ - Saa'J)niaanii3a' - J ^ alaa^jao-'al^o-'Oi/Ja (l) 

on a simple cubic lattice with two orbitals (labeled by a, (5) per site and intra-orbital nearest- 
neighbor hopping t. The Coulomb interaction matrix elements included here are the intra- 
orbital interaction [/, the inter-orbital density interaction V , and the inter-orbital exchange 
interaction (Hund's coupling) J. The last term represents the transverse part {S^'^S^p + 
S^S~p) of the Hund's coupling, and the density interaction term yields the longitudinal 
part Sl^Slp, so that altogether the Hund's coupling term has the form — JSj^.S^^. The 
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Hamiltonian therefore possesses continuous spin rotation symmetry, and hence the Goldstone 
mode must exist in the spontaneously broken symmetry state. 

Hund's couphng has been shown to strongly enhance ferromagnetism in an orbitally de- 
generate system by strongly suppressing the correlation-induced quantum corrections.^'' An 
effective quantum parameter [f/^ + {JV — 1) J^]/[f/ + (A/" — 1) J]^ was obtained for determin- 
ing the magnitude of quantum corrections for an A/'-orbital system, in analogy with 1/S 
for quantum spin systems. The rapid decrease of this quantum factor from 1 to l/JV as 
J/U increases from to 1 results in strong suppression of quantum corrections and hence 
significant stabilization of ferromagnetism by Hund's coupling, particularly for large J\f. 

In order to highlight the role of inter-orbital Coulomb interaction V and orbital fluctua- 
tions on magnetic couplings and excitations in this paper, we will first set J = 0. The case 
of finite Hund's coupling will be treated later in section VII. 

In a band ferromagnet, all information regarding carrier- induced spin interactions Jij = 
U'^4'ij and excitations are contained in the irreducible particle-hole propagator 0(q, u), which 
then yields the exact transverse spin fluctuation (magnon) propagator:— 



Our approach is to incorporate correlation effects in 0(q, u) in the form of self-energy and 
vertex corrections using a systematic expansion = 0^°) + 0^^) + 0^^) + ... which preserves the 
Goldstone mode order by order. Rooted within an inverse- degeneracy expansion scheme, this 
systematic approach is non-perturbative with respect to the interaction terms and therefore 
yields a controlled approximation which remains valid in the strong coupling limit. Con- 
tributions to the first order quantum correction 0*^^-' due to the Hubbard interaction U and 
Hund's coupling J have been discussed earlier.-Sfi 

The additional first order diagrams for arising from the inter-orbital interaction V are 
shown in Fig. 1. The diagrams shown here are for a saturated ferromagnet in which minority 
(D spin particle- hole fluctuations are absent. Here Fig. 1(a) and (b) represent quantum 
corrections due to electronic self energy renormalization by orbital fluctuations. Fig. 1(c) 
represents the corresponding vertex correction, and Fig. 1(d) represents vertex corrections 
involving coupling between transverse spin and orbital fluctuations. The vertex correction 
diagrams as in Fig. 1(d) (nine such diagrams) can be represented in terms of an effective 
spin-orbital interaction vertex Fgp.orb as shown in Fig. 2(a). The spin-orbital interaction 




(2) 




FIG. 1: First order diagrams for the irreducible particle- hole propagator arising from interaction 
V involving: (a,b) self energy corrections due to orbital fluctuations, (c) corresponding vertex cor- 
rection, and (d) vertex corrections involving spin-orbital coupling (altogether, nine such diagrams). 
The dashed and wavy lines represent interactions U and V ^ respectively. 




(a) 



— < ( < — 



(b) 



FIG. 2: The vertex correction diagrams in Fig. 1(d) can be represented (a) in terms of a spin-orbital 
interaction vertex Fsp-orb • The three diagrammatic contributions of Fsp-orb involving three- fermion 
vertices (b) generate the nine diagrams involving spin-orbital coupling. The missing fourth diagram 
vanishes because of the assumption of complete polarization. 
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vertex has three contributions involving three-fermion vertices, as shown in Fig. 2(b). The 
missing fourth diagram vanishes because of the assumption of complete polarization. 

As the Goldstone-mode condition C/0 = latg = Ois already exhausted by the zeroth- 
order (classical) term the sum of the higher order (quantum) terms (p^^^ + 0'-^-' + ... 
must exactly vanish at g = 0. For this cancellation to hold for arbitrary U, J, and V, 
each higher order term 0*^"^ in the expansion must individually vanish, implying that the 
Goldstone mode is preserved order by order. We will demonstrate this exact cancellation 
explicitly for the new contributions due to V in the first-order quantum correction 0*^^^ . 

Systematics in our two-orbital model can be formally introduced, in analogy with the 
inverse-degeneracy (l/A/") expansion for the Hubbard model, by: i) treating the two physical 
orbitals a, (3 as pseudo spins, ii) introducing M pseudo orbitals (/x) for each pseudo spin, 
and iii) generalizing the inter-orbital density interaction to {V/H) Ylii^iv''^ioiti'n'ipv Now, each 
interaction hue V yields a factor 1/J\f and each bubble yields a factor M from the summation 
over pseudo orbitals, resulting in an overall 1/M factor for the bubble series, and an overall 
(1/jV)" factor in the n*^-order quantum correction 0^"). 



III. SPIN AND ORBITAL FLUCTUATIONS 



The diagrammatic expansion above involves spin and orbital fluctuation propagators, 
the characteristic energy and momentum distribution of which are important in view of the 
spin-orbital coupling investigated in this work. The ladder series in Fig. 1(d) yields the 
effective intra-orbital transverse spin interaction: 

^-UxoKQM) l-f/xo(Q,^) Q + w^-^r/ 

in terms of the RPA-level magnon propagator Xsp-, having an advanced pole in the saturated 
ferromagnetic state {n^ — m, ^4. = 0). Here xo(Q,^^) is the bare antiparallel-spin particle- 
hole propagator, mq m and ujq are the magnon-mode amplitude and energy, and the 
small weight of gapped Stoner excitations has been neglected for simplicity. 

Similarly, the bubble series in Fig. 1 involving odd number of bubbles yields, in terms of 
the RPA-level orbital fluctuation propagator, the effective intra-orbital density interaction: 
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FIG. 3: Momentum-integrated orbital fluctuation spectral function (a) shows strong suppression of 
the para-orbiton energy scale near the orbital ordering instability; the low-energy integrated part 
(b) shows that low-energy orbital fluctuations are concentrated near (7r,7r,7r). 

near the orbital ordering instability where Vxot ~ 1- Here xot(Q)^) is the bare spin-'|" 
particle-hole propagator. The orbital fluctuation propagator is symmetric Xorb(— Q, — ^) = 
Xorb(Q,^), with a spectral representation: 

dn'imixoM^')] 



Xorb(Q,^^) = - / 

Jo 



(5) 



7r n — n' + ir] 

for its retarded part, with a continuum distribution over the orbital fluctuation energy Q'. 

Exactly at quarter fllling {m = n = 0.5 per orbital), the orbital fluctuation propagator 
diverges at Q = (vr, vr, vr) in the absence of any NNN hopping t' terms which destroy Fermi 
surface nesting, indicating instability towards staggered orbital ordering. With increasing 
hole doping, the spectral function peak shifts below (7r,7r,7r). 

Fig. 3(a) shows the momentum integrated orbital fluctuation (para-orbiton) spectral 
function ^Q(l/7r)ImXorb(Q; ^) with increasing interaction strength V. Here, and in the 
following, we have set the hopping term t = 1 = W/12 as the unit of energy, where W is 
the bandwidth. In analogy with the well-known para-magnon response with approaching 
magnetic instability, the para-orbiton energy scale is strongly suppressed from order band- 
width in the weak-coupling regime to relatively very low energies near the orbital ordering 
instability. 

Fig. 3(b) shows the momentum dependence of the low-energy {Q < 2) integrated part of 
the orbital fluctuation spectral function near the R point (7r,7r,7r), which shows that low- 
energy orbital fluctuations are concentrated near the wavevector (7r,7r,7r) corresponding to 
staggered orbital fluctuations. With increasing doping away from quarter fllling, the peak 
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shifts below (7r,7r, tt), indicating incommensurate fluctuations. 



IV. ORBITAL FLUCTUATION MAGNON SELF ENERGY 



We wiU flrst consider diagrams for in Fig. 1 (a,b,c) involving electronic self energy 
corrections due to orbital fluctuations and the corresponding vertex correction. Absent 
in the single-orbital case, these diagrams are characteristic of orbital degeneracy, inter- 
orbital interaction, and orbital fluctuations, and strongly influence magnetic couplings and 
excitations through electronic band renormalization, particularly in vicinity of the orbital 
ordering instability. The vertex correction diagrams (d) involving both spin and orbital 
fluctuations will be discussed in the next section. 

Quantum corrections to the irreducible particle-hole propagator in Eq. (2) yield the 
magnon self energy: 



in terms of which the magnon propagator x"+(q, a;) = m/[uj + coq — E(q, a;)]. The flrst- 
order magnon self energy corresponding to diagrams in Fig. l(a,b,c) involving only orbital 
fluctuations is then obtained by summing over the bosonic degrees of freedom of the orbital 
fluctuations: 



where r4 is the four-fermion vertex obtained by integrating out the fermionic degrees of 
freedom in the diagrams for shown in Fig. l(a,b,c). As will be further discussed in the 
following three subsections, this orbital fluctuation magnon self energy physically represents 
contributions due to (i) coupling between Stoner excitations and orbital fluctuations, and (ii) 
self-energy corrections involving band-energy renormalization and spectral-weight transfer. 

Using the spectral representation for the orbital fluctuation propagator (retarded part, 
since r4 has only advanced poles with respect to il), we obtain: 




(6) 




oo 



dQ' ImXorb(Q - q, ^0 
TT — uj — D,' + ir) 




(8) 
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where the average four-fermion vertex (r4)Q/ q/ denotes averaging over orbital fluctuation 
modes Q' = Q — q. Evaluation of the four-fermion vertex resolved into different contri- 
butions corresponding to distinct physical mechanisms, is discussed below. Term by term, 
the four-fermion vertex r4 vanishes identically for g = 0, in accord with the Goldstone mode. 



A. Stoner-orbital coupling 

In the single-orbital case, the magnon self energy due to spin-charge coupling included 
a Stoner contribution representing coupling of charge excitations with the gapped part of 
spin excitations.^^ In analogy, diagrams Fig. l(a,b,c) yield contributions which represent a 
Stoner-orbital coupling: 

rStoner^^2^/ I \ I I 1 



+n-trij V4-q-4 4-Q-4-Q+q + ^-^^. 

in which the first term represents the Stoner excitation mode (Q, Q) and the quadratic term 
is the Stoner-orbital interaction vertex, which involves only magnetic energy denominators. 
Here ~ ~ are the ferromagnetic state band energies for the two spins in terms of 
the free-particle energy = — 2t(cos + cos ky + cos k^) for the sc lattice and the exchange 
splitting 2 A = mil. The band energy superscripts +(— ) refer to particle (hole) states. 
There is no restriction on the energy ej^^Q^^ in Eq. (9) as both particle and hole states 
contribute. Of the four terms in this quadratic interaction vertex, the two square terms 
arise from diagrams Fig. 1(a) and (b), while the cross terms are from diagram Fig. 1(c); the 
characteristic quadratic structure therefore stems from orbital fluctuations renormalizing 
electrons of both spins, and is clearly absent in the single-orbital case involving spin-down 
renormalization only. This intrinsic quadratic structure resulting from orbital degeneracy 
yields a characteristic non-Heisenberg (1 — cosg)^ magnon self energy, resulting in no spin 
stiffness reduction but strong zone-boundary magnon energy reduction. 



B. Electronic band renormalization 



Due to exchange of inter-orbital fluctuations in the diagram Fig. 1(a) involving intermedi- 
ate spin- 'I' states, the spin-'|' hole (particle) energies are pulled down (pushed up), increasing 
the particle-hole energy gap, and thereby suppressing the particle-hole propagator (p. In- 
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eluding the eorresponding eontributions from the vertex correction diagram Fig. 1(c), we 
obtain the electronic band renormalization contribution: 

rr = -u^J: U ^3^^) ( ^.. - J. ^ ) ' (10) 

involving one charge and two magnetic energy denominators. The finite infinitesimal term 
ir) as in Eq. (9) has been dropped for compactness. 

C. Spectral weight transfer 

The electronic self energy correction in diagram Fig. 1(a) also results in spectral- weight 
transfer and redistribution between occupied and unoccupied spin-t states. However, there is 
no net change in occupancy and magnetization. The corresponding spectral weight transfer 
contribution: 

pspectral ^ -U^ i ^fT ^t" , Q ~ ) ( M J^TT, ~ M TTT ) ^^^^ 

k Vek-Q+q-^k \ek_q-ek +^ ^k-Q - ^k-Q+q + 

involves one magnetic and two charge energy denominators. The first (negative) contribution 
corresponds to loss of spin-t hole spectral weight due to transfer to particle states, and the 
second (positive) contribution corresponds to the reverse process. 

Cancellation of most singular contributions 

Singular contributions in Eqs. (9-11) for r4 exactly cancel out. For example, the most 
singular contribution in Eq. (9) involving two powers of the vanishing energy denominator 
{ej^Q — e]^'^Q_|.q -I- uj) exactly cancels with the corresponding contribution from Eq. (10). 
Similarly, the next most singular contributions in Eqs. (9-11) also exactly cancel out. 

Average over orbital fluctuation modes 

The average (r4)Q/^f2/ of the four-fermion vertex over orbital fluctuation modes directly 
yields the magnon self energy from Eq. (8). Since orbital fluctuations peak below (tt, 7r,7r) 
for finite doping [Fig. 3], the vertex (r4) was estimated by averaging over a selected Q' 
region (—0.8 < cosQ'^ < —0.6) near (7r,7r, tt), with ImXorb(Q', ^0 assumed fiat inside and 
zero outside. 
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FIG. 4: The Stoner, band, and spectral contributions of the averaged four-fermion vertex (r4) 
shows strong non-Heisenberg (1 — cos (7)^ behaviour in the T-X direction. 

The q dependence of the averaged four-fermion vertex (r4(g)) in the F-X direction of the 
Brillouin zone is shown in Fig. 4. The band contribution is negative due to renormahzation 
of band energies by orbital fluctuations, as discussed above. The band contribution rapidly 
diminishes in the strong coupling limit [Fig. 4(b)], as does the Stoner contribution. The 
spectral contribution survives in the strong coupling limit. All contributions have strongly 
non-Heisenberg character, with negligible magnitude at small q but rising sharply near 
the zone boundary, implying no spin stiffness correction but appreciable zone boundary 
magnon softening. Orbital fluctuation modes near (tt, 0, tt) etc. also yield significant non- 
Heisenberg character to the interaction vertex F4(g), and the three contributions exhibit 
similar behaviour. 

How does the orbital fluctuation magnon self energy compare with the bare magnon 
energy? Taking the orbital fluctuation energy Vt' to be negligible in comparison to the 
bandwidth near the orbital ordering instability, and the estimated average (r4(q))Q/ 0.3 
near the zone boundary from Fig. 4(c), we obtain (for m = 0.35 and V = 3): 

Sorb(q) ~ m(\/V2)(F4(q))Q' = 0.35 x (9/2) x 0.3 ^ 0.5 (12) 

which is comparable to the bare zone-boundary magnon energy ^ 0.35(1 — cosg) ~ 0.7 
for realistic strength of the inter-orbital interaction. 

V. SPIN-ORBITAL COUPLING MAGNON SELF ENERGY 

In the previous section, we considered the diagrams of Fig. 1 (a,b,c) involving electronic 
self-energy corrections due to orbital fluctuations. We will now examine the vertex correction 
diagrams of Fig. 1(d) representing spin-orbital coupling, which are particularly important 
near the orbital ordering instability where orbital fluctuations are soft. In contrast to the 
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single-orbital case where self energy and vertex correction diagrams were of qualitatively 
similar order,— introduction of the orbital degree of freedom lifts this constraint and allows 
qualitatively independent self energy and vertex correction contributions. The corresponding 
first-order magnon self energy: 



''sp— orb 

= — m 



Q 

^[xsp(Q,fi)][rsp-orb(Q,q,fi,w)]^[xorb(q- Q,w - ^)] (13) 

where Xsp and Xorb are the spin and orbital fluctuation propagators (section III), and 
Tsp-orb = U^V^3 represents the spin-orbital interaction vertex in terms of the three-fermion 
vertex Fs, evaluation of which is discussed in the Appendix. 

In analogy with the spin-charge coupling process,— this correlation-induced spin-orbital 
coupling is analogous to a second-order Raman scattering process in which the magnon (q, u) 
scatters into an intermediate-state magnon (Q, Q) along with an internal orbital excitation 
(q — Q, — f2), leading to significant magnon energy renormalization and magnon damping. 

The Q, Q integrals in the above equation represent integrating out the bosonic (both spin 
and orbital) degrees of freedom. As the magnon propagator [xsp(Q, ^)] is purely advanced 
in nature, only the retarded part of the product [rsp_orb]^[Xorb] contributes in the Q integral. 
For simplicity, considering the dominant contribution to the spectral representation of this 
product coming from the orbital fluctuation propagator, from Eqs. (13), (3), (5), and the 
symmetry property given above Eq. (5), we obtain: 



^lD-orb(q'^) = 



sp-orbVH,-; - -n.^ J — y^^^o^_^^^ [r.p-orb(Q,q,f^,c.)]2^ ^ n-u-Q' + zr^ 
m^ir'-^ W.b(Q - q, n') (14) 



+ oo + \L' — IT] 

Due to magnon decay into internal spin and orbital excitations (— = Wq = Uq + Q'), the 
above magnon self energy yields a finite imaginary part representing finite magnon damping 
and linewidth, as discussed earlier for spin-charge coupling.— 

An approximate evaluation of the resulting spin-orbital magnon self energy illustrates the 
importance of the orbital fluctuation softening near the orbital-ordering instability. With 
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Qorb — ^' and ilgpin = representing characteristic orbital and spin fluctuation energy 
scales, we obtain (for u — 0): 

E^^L (q) ~ ;„2 ([rsp-orb(q)]^)Q',n' ^^g^ 

where the angular brackets ( ) again refer to averaging over the orbital fluctuation modes 
Q' = Q — q, as in Eq. (8). Far from the orbital ordering instability, the orbital fluctuation 
energy Qorb is of order bandwidth W, which strongly suppresses the magnon self energy. 
However, near the orbital-ordering instability, spin-orbital coupling becomes important due 
to the relatively much smaller energy denominator Jlgpin + ^^orb ~ t- 

The spin-orbital interaction vertex rsp_orb is obtained by integrating out the fermion 
degrees of freedom in the three-fermion interaction vertices. This interaction vertex explicitly 
vanishes at momentum g = in accordance with the Goldstone mode requirement, and 
yields the dominant q dependence of the magnon self energy. In order to illustrate the 
characteristic non-Heisenberg character of the interaction vertex, we consider its magnetic 
part with energy denominators involving the Stoner gap. This term qualitatively differs from 
the charge part of the vertex with energy denominators involving excitation energies of order 
bandwidth. Evaluation of the three-fermion vertices contributing to Egp-orb is discussed in 
the Appendix. 

For the purely magnetic part, we obtain: 



4-+ t- I ■ 1+ t I 



where there is no restriction on the e]^_,_q_Q energies. The leading order q dependence of 
the above term is approximately (1 — cosq) as shown below. For a fixed orbital fluctua- 
tion momentum Q' = Q — q, with Q' close to (tt, tt, tt) corresponding to staggered orbital 
fluctuations, expansion in powers of e/A of the square-bracket terms in Eq. (16) yields a 
contribution (ek_q — ek)/2A ~ (1 — cos q) to leading order. The formally second-order struc- 
ture of the spin-orbital coupling magnon self energy (Eq. 15) involving [Esp.orb]^ therefore 
directly yields an intrinsically non-Heisenberg (1 — cosg)^ contribution to the magnon self 
energy, which yields no spin stiffness quantum correction, but significant magnon energy 
reduction and damping near the zone boundary, and therefore accounts for the observed 
zone-boundary magnon anomalies. 
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FIG. 5: Momentum dependence of (a) the magnetic part [r™^^jj.jj]^ and (b) the total spin-orbital 
interaction term [Fsp-orb]^ shows a strong anomalous behaviour in the F-X direction. While stag- 
gered orbital fluctuation modes near (vr, tt, vr) yield the dominant contribution to the magnetic part 
(a), it is modes near (7r/2, 7r/2, 0) which yield the dominant contribution to the total (b), which is 
also enhanced with hole doping as in (a). 

The q dependence of the magnetic part [r^'^^j.j^]^ of the spin-orbital interaction term is 
shown in Fig. 5(a). Being small in comparison to the Stoner gap, the boson energies u and 
Q in Eq. (16) were set to zero, and the momentum Q' was selected in a range near (tt, tt, tt), 
as in Fig. 4. While [r^'^^j.^,]^ shows strong anomalous momentum dependence for orbital 
fluctuation modes near (vr, vr, vr), when both magnetic and charge terms (see Appendix) are 
included, the net contribution from these modes is small due to a cancellation, as shown in 
Fig. 5(b). 



VI. ORBITAL FLUCTUATIONS NEAR (7r/2, 7r/2, 0) 

In contrast, for orbital fluctuation modes with momentum Q' near (7r/2, 7r/2, 0), the 
above cancellation is avoided as the magnetic contribution to Fsp_orb is small. We flnd that 
the total spin-orbital interaction term [Fsp_orb]^ exhibits a strong anomalous momentum 
dependence [Fig. 5(b)], and that it is strongly enhanced with increasing hole doping, as the 
band filling changes from 0.5 (quarter filling) to 0.25 (one-eighth filling). 

In order to quantitatively examine the effect of this anomalous momentum behaviour 
of rgp_^jj.j^ on the magnon dispersion, the magnon self energy was evaluated approximately 
using Eq. (15). On averaging [Fsp_orb]^ over Q', we find that the anomalous momentum 
behaviour remains qualitatively similar for \Q'^\ ^ 1 and drops sharply for Q'^^ 1, whereas in 
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FIG. 6: The renormalized magnon energy due to spin-orbital coupling [Eq. (17)] shows strong zone- 
boundary magnon softening in the T-X direction, which becomes more pronounced with increasing 
hole doping. 

the Q'^ — Q'y plane the dominant and qualitatively similar contribution comes from diagonal 
strips along + Q^l = tt, yielding a factor of ~ 1/2 on planar averaging, resulting in an 
overall phase-space factor of ~ (l/3)(l/2). 

The magnon self energy was therefore estimated using Eq. (15) with rgp_Qj.jj as obtained 
in Fig. 5(b) for the orbital fluctuation mode (7r/2, 7r/2, 0), and including the above phase- 
space factor (taken as 1/10 in our calculations) to account for momentum averaging over 
orbital fluctuation modes. Also, the bare magnon energy cjq and the orbital fluctuation 
energy f2orb were taken to be of order t. The renormalized magnon energy: 



is shown in Fig. 6 for different band fillings n = (1 — x)/2. It is evident that while the 
spin stiffness remains essentially unchanged, the magnon self energy at the zone-boundary 
becomes comparable to the bare magnon energy for realistic values of the inter-orbital in- 
teraction V ^ resulting in a zone-boundary softening which is strongly enhanced with hole 
doping. In the presence of staggered orbital correlations such as near (7r/2, 7r/2, 0), weighted 
averaging over Q' with a peaked orbital spectral function will further enhance ([Tsp^orb]^) 
and therefore the anomalous magnon self energy. 



,0 v(i) 



orb(q) 



(17) 



17 




FIG. 7: Period 4a planar staggered orbital correlations corresponding to orbital fluctuation modes 
with momentum near (7r/2, 7r/2, 0). 



The orbital fluctuation modes near (7r/2, 7r/2, 0) correspond to period 4a planar staggered 
correlations, as shown in Fig. 7. Such orbital correlations have been observed in the CE-type 
charge-ordered phase of the half-doped (x = 0.5) manganites such as the narrow-band com- 
pounds like Pri_xCaxMn03 and Lai_xCaxMn03,— and the layered material Lai/2Sr3/2Mn04 
in which magnetic excitations were found to be dominated by ferromagnetic couplings.— 
Therefore, as x approaches 0.5, orbital fluctuation modes near (7r/2, 7r/2, 0) may form the 
dominant contribution to the low-energy part of the orbital fluctuation spectral function. 

The different layers in Lai/2Sr3/2Mn04 are magnetically decoupled due to negligible in- 
terlayer couplings. Extension of the present investigation of spin-orbital coupling to the 
two-dimensional case is therefore of interest with respect to renormalization of magnetic cou- 
plings. In the case of spin-charge coupling, the renormalized magnon dispersion for a square 
lattice does show strong zone-boundary softening near (vr, 0) and (0,7r), while the magnon 
energy near (vr, vr) remains undiminished, indicating softening of the nearest-neighbour fer- 
romagnetic bonds but strong ferromagnetic correlations along the diagonal directions.— 

Orbital fluctuation modes near Q' = (ivr, 0, 0) and (0, ±7r, 0) were also found to yield 
signiflcant contribution to the spin-orbital interaction vertex. Involving similar density (njQ,+ 
riip) on all sites, such conflgurations should, however, be relatively suppressed by an intersite 
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density interaction VniTij. In contrast, allowing for reduced density on alternating "empty" 
sites with vanishing orbital "magnetization", configurations corresponding to modes near 
(7r/2, 7r/2, 0) with period 4a orbital correlations minimize the intersite interaction energy 
VniTij, and would therefore be relatively more important. 



VII. FINITE HUND'S COUPLING 



So far we had set the Hund's coupling J = in order to highlight the role of inter-orbital 
interaction and fiuctuation. For finite J, it is convenient to proceed in two steps. The part 
^ = J of the inter-orbital interaction V together with Hund's coupling effectively amounts to 
a purely magnetic interaction — JSj^.S^/?, and has been investigated earlier.— This is because 
for V = J, the inter-orbital interaction (J — J5cra')niaanii3a' acts only between opposite-spin 
electrons, and so the resulting diagrammatics is similar to the Hubbard interaction case. 
The remaining part {V — J) is purely non-magnetic, and yields diagrammatic contributions 
essentially as in section II. Thus, corresponding to Fig. 1 (a,b,c) diagrams involving the 
bubble series for the orbital fiuctuation propagator, the magnon self energy is obtained by 
simply replacing V hj (V — J) in Eq. (7). 

For the spin-orbital coupling magnon energy, however, the interaction ladders in diagrams 
Fig. 1(d) and Fig. 2 now involve Hund's coupling J as well, and since the transverse part 
of J fiips the orbital index in the ladder series, there are now two contributions to the 
irreducible particle-hole propagator: 

0„^(q,u;) = -V / ^U:i{Cl,Q)[r,{Cl,ci,Q,i^)?V:/{q-Q,u-Q) (18) 

involving intra (yU = a) and inter {fi = /3) orbital spin and orbital fiuctuations. Involving 
ladders of U and J, the effective intra and inter-orbital transverse spin interactions: 



u+ ^ u- 



Xo _^ Xo 
1 - U+xo 1 - U-xo 



2 . .... rr . (19) 

can be expressed as in-phase (/i = a) and out-of-phase (/i = (3) combinations of the acoustic 
and optical branches,— with U"^ = U ± J. Similarly, the effective intra and inter-orbital 
density interactions: 

{V-Jfxot 



^eff 



1 - - jyxi-^ 



~ i-{v- JYxl, ^ ^ 
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involve odd and even number of bubbles. 

Now, in the investigation of role of Hund's coupling on quantum corrections,^*^ it was 
shown that the inter-orbital component t/^(Q,f2) yields a small ^-integrated contribution 
due to partial cancellation from the out-of-phase combination of the acoustic and optical 
modes, essentially reflecting an inter-orbital incoherence. Therefore, the inter-orbital com- 
ponent (pctji should be relatively small, and the intra-orbital component ipaa is approximately 
given by Eq. (15), with V replaced by {V — J) and U replaced by {U + J). 

Due to the purely opposite-spin density interaction (J — JSo-a')nia(Tni/3a', the interaction 
line U connected to the bubble in the third diagram in Fig. 2(b) for the three-fermion vertex 
Fa is also replaced by {U + J). This is consistent with the enhanced exchange splitting to 
{U + J)m in the energy denominator, which changes the magnon pole condition to 
{U + J)x^ = 1; and ensures that the three-fermion vertex F3 exactly vanishes at g = in 
accordance with the Goldstone-mode condition (see Appendix). 

VIII. EXTENSION TO THE FERROMAGNETIC KONDO LATTICE MODEL 

Magnetic couplings and excitations in ferromagnetic manganites have been theoretically 
investigated using the ferromagnetic Kondo lattice model (FKLM) and its strong coupling 
limit, the double exchange ferromagnet. In this model, the S = 3/2 core spins due to 
localized t2g electrons of the magnetic Mn"*""*" ions are exchange coupled to the mobile eg 
band electrons, represented by an interaction term — J' Sj.cTj, with J' of order bandwidth 
W. The fermionic representation approach for evaluating magnon self energy corrections 
in the FKLM, which allows conventional diagrammatic tools to be employed for evaluating 
quantum corrections beyond the leading order,'^ can be readily extended to include effects 
of orbital fluctuations. 

With an inter-orbital interaction V included between the two degenerate eg orbitals, the 
FKLM magnon self energy due to spin-orbital coupling can be directly obtained from our 
magnon renormalization analysis of Fig. 6. The required correspondence is: U ^ J and 
m — 7- 25*, so that the exchange splitting Um ^ 2JS. The FKLM magnon self energy 
is thus obtained from our calculated Hubbard model result using a multiplicative factor 
/ = J^{2Sf/U^m^. With U/t = 20, m = 0.35, J /t = 4, and 2S = 3, we obtain / ^ 1/6. 
As the FKLM bare magnon energy ~ (t/18)(l — cosg) is smaller than the Hubbard model 
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bare magnon energy ~ (t/3)(l — cosg) by roughly the same factor, the FKLM renormahzed 
magnon energy = 0;° — Ssp_orb(ci) will also be as in Fig. 6, only scaled down by a factor 
(1/6). Taking the hopping energy scale t ~ 200meV corresponding to a realistic bandwidth 
~ 2eV, the magnon energy scale in Fig. 6 is ~ 30meV for ferromagnetic manganites, in 
agreement with the measured magnon energies. -- 

Since the FKLM magnon self energy goes as the fourth power of J' /t explicitly, the 
anomalous zone-boundary softening effect should be especially pronounced in narrow-band 
systems. Indeed, zone-boundary magnon softening is clearly seen to occur in the relatively 
low-Tc or narrow-band materials, and broad-band materials such as Lao.yPbo.sMnOs show 
nearly Heisenberg behavior, in agreement with this prediction. Systematic studies of doping 
dependence of the zone-boundary magnon softening indicates that the higher the doping 
level, the larger the zone-boundary softening.— Furthermore, doping dependence of spin 
dynamics indicates that the measure spin stiffness D ~ 160 ± 15 meV remains essentially 
unchanged, while the zone-boundary magnon softening (denoted by the ratio J4/ Ji) is found 
to be enhanced linearly with increasing doping. 

IX. CONCLUSIONS 

The correlated motion of electrons in the presence of strong orbital fluctuations and 
correlations was investigated in an orbitally degenerate ferromagnet with two orbitals per 
site with respect to magnetic couplings and excitations. A systematic Goldstone-mode- 
preserving approach was employed to incorporate correlation effects in the form of self 
energy and vertex corrections, so that both long-wavelength and zone-boundary magnon 
modes could be studied on an equal footing. Our investigation focussed on the anomalous 
momentum dependence of the three- and four-fermion interaction vertices which determine 
the magnon self energy and the role of different orbital fluctuation modes, particularly near 
the orbital ordering instability where orbital fluctuations are relatively soft. 

Orbital fluctuations were generically shown to impart an intrinsically non-Heisenberg 
(1 — cosg)^ character to the magnon self energy in the F-X direction of interest. This 
generic behaviour was ascribed to a quadratic structure of the spin-orbital interaction vertex 
resulting from a new class of diagrammatic contributions associated with the orbital degree 
of freedom. These diagrams are absent in the single-orbital case, the essential difference 
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being that orbital fluctuations couple to electrons of both spin. 

The absence of any contribution in this non-Heisenberg magnon self energy for small 
q implies that the spin stiffness is not renormalized by orbital fluctuations generated by 
the inter-orbital density interaction V . In a multi-orbital ferromagnet, with intra-orbital 
Coulomb interaction [/, inter-orbital interaction V , Hund's coupling J, and orbital degener- 
acy A^, the spin stiffness therefore continues to be essentially determined by the intra-atomic 
factors f/, J, A/", through the effective quantum parameter [f/^ -|- (A/" — 1) J^]/[f/ + (A/" — 1) J]^ 
as obtained earlier,— and the interaction V does not play an important role. 

However, the strong enhancement of magnon self energy near the zone boundary resulted 
in a strong anomalous magnon softening in the F-X direction, which increases significantly 
with hole doping away from quarter filling. Our investigation thus clarifies the completely 
different roles of interactions J and V , representing magnetic and charge parts of the inter- 
orbital Coulomb interaction. While Hund's coupling enhances ferromagnetism by strongly 
suppressing the effective quantum parameter, orbital fluctuations and correlations due to 
V destabilize ferromagnetism by strongly suppressing zone-boundary magnon energies near 
the orbital ordering instability. 

With regard to relative importance of different orbital fluctuation modes, staggered fluc- 
tuations with Q near (vr, vr, vr) and (vr, 0, vr) were found to be most important for the orbital 
fluctuation magnon self energy. The spin-orbital coupling magnon self energy was found to 
be strongly sensitive to orbital fluctuation modes due to an interference between magnetic 
and charge terms in the interaction vertex. Thus, while the magnetic part showed strong 
anomalous momentum dependence for orbital modes near (7r,7r, vr), the net contribution to 
the total vertex was found to be small due to cancellation with the charge term. Rather, 
fluctuation modes near (vr/2, 7r/2, 0) were found to be important for the total vertex includ- 
ing the charge part. The strong zone-boundary magnon softening near (0, 0, vr), arising from 
staggered orbital fluctuations with Q near (7r/2, 7r/2, 0), suggests an instability towards a 
composite structure of spin-orbital correlations involving period 4a orbital ordering in fer- 
romagnetic planes and intra-orbital AF spin correlations in the perpendicular direction. 

These results provide a plausible explanation of the observed anomalies in neutron scatter- 
ing studies of spin-wave excitations in ferromagnetic manganites, where spin stiffness is seen 
to remain essentially unchanged whereas the zone-boundary magnon softening is enhanced 
with increasing hole doping and the approach towards CE-type charge-orbital ordered states 



22 



near x = 0.5. Our results of strong anomalous magnon self energy contribution from dif- 
ferent orbital fluctuation modes such as (7r,7r,7r) and (7r/2, 7r/2, 0) show the zone-boundary 
softening to be a more generic feature of spin-orbital coupling. Only ferromagnetic orbital 
correlations extremely close to the orbital ordering instability were found to yield significant 
magnon softening in earlier studies.— 

The observed zone-boundary magnon softening has been usually modeled by including a 
fourth neighbour interaction term J4.^- As J4 yields no contribution to the zone-boundary 
magnon energy, but contributes significantly to spin stiffness, it must be accompanied by a 
corresponding reduction AJi = 4J4 in the NN coupling so that the spin stiffness remains 
unchanged, as observed experimentally; the net magnon energy reduction then has the non- 
Heisenberg form 2 74(1 — cos g)^. Our anomalous magnon self energy result of this form thus 
provides fundamental insight into the role of orbital fluctuations on magnetic couplings and 
excitations. 

Instead of the inter-orbital interaction VniaUi/j considered here, inter-site interactions 
VniUj would generate similar diagrammatic contributions to the magnon self energy, which 
become important near the charge-ordering instability where charge excitations become 
relatively soft, resulting in similar magnon self energy and anomalous zone-boundary magnon 
softening. 
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Appendix 



The spin-orbital interaction vertex is obtained as: 



sp— orb 



(21) 



in terms of the three types of fermion vertices shown in Fig. 1(e), which are evaluated by 
integrating out the fermion degrees of freedom as discussed below. 

For the first fermion vertex with interaction line V attaching to spin-J, fermion, we obtain: 
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and a magnetic-charge term: 
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(24) 

=k-Q -^k-Q+q + ^-^^y V^k-Q+q-^k -UJ + n-tr]^ 

In Eq. (23), there is no restriction on the fermion energy ej^„Q_,_q as contributions with both 
particle (+) and hole (-) energies are included. Finally, the third vertex Fg'^'' involves a U 
interaction line and a spin-']" bubble attached to the spin-J, fermion lines, and is given by: 



.3 --ri^'^ uxoM-Q,^-^) (25) 

For q = 0, the spin-orbital interaction vertex Fsp^orb vanishes identically, ensuring that the 
Goldstone mode is explicitly preserved. Both —T^^^ and vf^^ reduce to xo(Q)^)/(2A + oo) 



for q = 0, whereas both —Tf^^ and T''^' reduce to Xot(Q>'^ ~ ^) setting UxoiQ,^) = 1 
at the magnon pole, so that from Eq. (21) Fsp-orb = for g = 0. 

Singular contributions in the three- fermion vertex F3 exactly cancel out, as for the four- 
fermion vertex F4 discussed below Eq. (11). Thus, the singular contribution of Eq. (23) 
due to vanishing energy denominator e^-Q ~ ^k-Q+q + exactly cancels the corresponding 
contribution from the second term of Eq. (24). 
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